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Abstract 

In a recent paper [Phys. Rev. Lett. 93, 130401 (2004)], we proposed the idea of expanding 
the space of variations in variational calculations of the energy by considering the approximate 
wave function ip to be a functional of functions x '■ ip = V4x] rather than a function. The space of 
variations is expanded because a search over the functions x can i n principle lead to the true wave 
function. As the space of such variations is large, we proposed the constrained-search — variational 
method whereby a constrained search is first performed over all functions x such that the wave 
function functional ?/>[x] satisfies a physical constraint such as normalization or the Fermi-Coulomb 
hole sum rule, or leads to the known value of an observable such as the diamagnetic susceptibility, 
nuclear magnetic constant or Fermi contact term. A rigorous upper bound to the energy is then 
obtained by application of the variational principle. A key attribute of the method is that the wave 
function functional is accurate throughout space, in contrast to the standard variational method 
for which the wave function is accurate only in those regions of space contributing principally 
to the energy. In this paper we generalize the equations of the method to the determination of 
arbitrary Hermitian single-particle operators as applied to two-electron atomic and ionic systems. 
The description is general and applicable to both ground and excited states. A discussion on 
excited states in conjunction with the theorem of Theophilou is provided. Here we construct 
new analytical 3-parameter ground state wave function functionals for the negative ion of atomic 
Hydrogen and the Helium atom through the constraint of normalization. We present the results 
for the total energy E, the expectations of the Hermitian single-particle operators W = rf,n = 
—2, —1, 1, 2, W = Yli <5( r i)) an d W = J2i <5( r i— r )> the structure of the nonlocal Coulomb hole charge 
p c ( rr ') as a function of electron position r, and the expectations of the two particle operators 
u 2 , u, 1/u, 1/u 2 , where u = |r$ — Tj\. The results for all the expectation values are remarkably 
accurate when compared with the 1078-parameter wave function of Pekeris, thereby demonstrating 
the accuracy of the wave function functionals throughout space. The same accuracy is exhibited 
for the Coulomb holes when compared with the 'exact' ones determined by Slamet and Sahni. The 
point of the improved spatial accuracy of these functionals is further made by comparisons with the 
results of other wave functions that are not functionals such as those due to Hartree-Fock theory 
and Caratzoulas-Knowles. We conclude by describing our current work on how the constrained- 
search — variational method in conjunction with Quantal density functional theory is being applied 
to many-electron atoms. 
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I. INTRODUCTION 

In recent work jl| , we proposed the idea of expanding the space of variations in standard 
variational calculations of the energy Q] , thereby allowing for an improvement of the energy 
in such calculations. Equivalently, a required level of accuracy could be achieved with 
fewer variational parameters. In the traditional application of the variational principle, the 
space of variations is limited by the choice of analytical form for the approximate wave 
function. For example, if Gaussian or Slater-type orbitals or a linear combination of such 
orbitals is employed in the energy functional, the variational space is limited by this choice 
of functions. The proposed manner by which the space of variations can be expanded is 
by considering the wave function ^ to be a functional of a set of functions \ : ip = ip[x], 
rather than a function. This permits a greater flexibility for the wave function ip[x] because 
the functions x rnay be chosen such that ip[x] reproduces any well-behaved function. In 
principle, a search over such functions can lead to that function x for which i/j[x\ is the true 
wave function. 

The space over which the search for the functions x is to be performed, however, is 
simply too large for practical purposes, and a subset of this space must be considered. We 
define the subspace over which the search for the functions x is to be performed by the 
requirement that the wave function functional ip[x] satisfy a constraint. Typical constraints 
on the functional ip[x] are those of normalization, the satisfaction of the Fermi-Coulomb 
hole sum rule, the requirement that it lead to observables such as the electron density, 
nuclear magnetic constant, diamagnetic susceptibility, Fermi contact term, or any other 
physical property of interest. With the wave function functional ip[x] thus determined, a 
rigorous upper bound to the energy is obtained by application of the variational principle. 
In this way, not only is a particular property of interest or constarint obtained exactly, the 
energy is also determined accurately since the variational principle ensures it is correct to 
second order in the accuracy of the wave function. We refer to this method of determining 
an approximate wave function as the constrained- search — variational method. The method 
is general in that it is applicable to both ground and excited states. 

An attribute of constructing a wave function functional ip[x] via the constrained-search — 
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variational method is that there is an improvement in the structure of the wave function 
throughout all space. Thus, both single-particle expectations representative of different 
parts of space as well as two-particle expectations involving two different points in space 
are obtained accurately. 

As in standard variational calculations, the satisfaction of constraints imposed on the 
wave function functional, while ensuring the exactness of a specific property or properties, 
will nonetheless lead to a less accurate upper bound to the energy provided the space 
of variations remains fixed. Any such decrease in the accuracy of the upper bound can, 
however, be offset by an increase in the space of variations. 

The concept of the wave function ^ as a functional ip[x] is general in that the space of 
variations may be expanded through the functions x- The number of functions x are al so 
independent of the electron number N. This contrasts with the Hartree-Fock theory [^] 
Slater determinant $[</>i] wave function which is also a functional but one of the iV-electron 
spin-orbitals <pi. Furthermore, there is no variational-flexibility of these spin-orbitals once 
they have been determined self-consistently by solution |4| of the Hartree-Fock equations. 
The space of variations cannot be expanded further, and therefore the Hartree-Fock theory 
wave function functional cannot be adjusted via the spin-orbitals 0, to be the true wave 
function. Thus, this wave function functional constitutes a point in the variational space as 
defined for the functional ip[x]- The determinantal functional $[0j] is therefore not general 
in the manner of the proposed ip[x]- 

In our original work [1] we had noted that the constrained-search — variational method 
could be extended to the determination of arbitrary Hermitian single-particle opera- 
tors. In sect. 2 we present the equations of this generalization as applied to the ground 
state of the negative ion of atomic Hydrogen, the Helium atom, and its isoelectronic 
sequence. The extension of these ideas to excited states in conjunction with the theorem 
of Theophilou P| is also described. We had also indicated various ways by which the 
results presented in our prior work could be improved. One such mechanism was to 
improve the prefactor in the correlated-determinantal wave function functional. In sect. 3 
we present the results of the application of the method with such an improved 3-parameter 
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analytical wave function functional to the ground state of the negative ion of atomic 
Hydrogen and the Helium atom, with normalization as the constraint. We present the 
results for the total energy E, the expectations of the Hermitian single-particle operators 

W = J2i r i, n = -2,-1,1,2, W = J2i 5 ( r i)i and W = Ei 5 ( r i - r )> the structure of the 
dynamic (nonlocal) Coulomb hole charge p c (rr') as a function of electron position r, and 
the expectations of the two particle operators u 2 ,u,l/u,l/u 2 , where u = |r$ — Tj\. The 
results for all the expectation values are remarkably accurate when compared with the 
1078-parameter wave function of Pekeris p, thereby indicating the accuracy of the wave 
function functionals throughout space. The same accuracy is exhibited in a different way by 
the comparison of the Coulomb holes with those of the essentially exact holes determined 
by Slamet and Sahnij?]]. The results for the energy and two particle expectations are 
far superior to those of Hartree-Fock theory as expected. However, the single-particle 
expectations are essentially equivalent since such expectations within Hartree-Fock theory 
are correct to second order The comparison with Hartree-Fock theory demonstrates 
how two square-integrable normalized antisymmetric wave functions can lead to essentially 
the same electron density^, but that one can be significantly superior to the other. Our 
results are also superior to those of the 3-parameter variational Caratzoulas-Knowles wave 
function that has a similar correlation term as ours but is not a functional. In the 
concluding section 4, we describe our current work on how the ideas of constructing wave 
function functionals are being applied in conjunction with Quantal density functional 
theory to the many-electron atom. 



II. CONSTRAINED- SEARCH VARIATIONAL METHOD 

In this section we present the generalization of the constrained-search — variational 
method for constraints whereby typical observables such as the diamagnetic susceptibil- 
ity, nuclear magnetic constant, Fermi contact term, and the constraint of normalization 
are determined exactly. For the two-electron systems represented by the negative ion of 
atomic Hydrogen, the Helium atom, and its isoelectronic sequence, these properties are rep- 
resented by the expectations of the single-particle operators W = r\ + r|, W — l/r 1 + l/r 2 , 
W = 5(ti) + 5(ra), and W — 1. For these two-electron systems, the Hamiltonian in atomic 
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units (e = h = m = 1) 

H = -\v\-bl- --- + -, (1) 
2 2 ri r 2 r i2 

where ri, r 2 are the coordinates of the two electrons, ri 2 is the distance between them, and 
Z is the atomic number. We next choose the form of the wave function functional to be of 
the general form 

il>[ X } = Hs,t,u)[l-f( X ;s,t,u)}, (2) 
with t,u) a pre-factor and f(x'i s,t,u) a correlated correction term: 

f(s, t, u) = e-^l + qu) [1 - X (g; s, t, u) (1 + u/2)], (3) 



12j, and q is a 



where s = r% + r 2 , t = r± — r 2 , w = ri 2 , are the Hylleraas coordinates 
variational parameter. Note that any two-electron wave function in a ground or excited 
state maybe expressed in this form. The key to the wave function functional is the 
determination of the functions x{q\ s,t,u). The prefactor may be chosen to be of some 
analytical form with variational parameters as in the present work, or the Hartree-Fock 
theory wave function jj], or determined self-consistently within the framework of Quantal 
Density Functional Theory 



For purposes of clarity, and thereby of subsequent analytical ease of solution, we assume 
the prefactor to depend only on the variables s and t: $ = $(s,t), and for the ground l 1 ^ 
state to be of the analytical form 

s,t) = Ne- as cosh{(3t) = ^.[ e - z ^ e ~ Z2r2 + e~ z ^e~ z ^\ (4) 

where different orbitals are allocated to electrons with up and down spins, a and (3 are 
variational parameters, Z\ — (a — 0) , Z 2 = (a + (3) , and iV is the normalization constant 
(See the Appendix). (Note that the normalization of the prefactor is independent of 
that of the wave function.) We further assume that x is a function only of the variable 
s: \P = ty[x(q;s)}. (The space of variations could be expanded further by assuming the 
function x to depend additionally upon the variable t, or still further by a dependence on t 
and u as well.) 



6 



The wave function functional ^[x{q] s)] for the ground state then satisfies the electron- 
electron cusp condition which in integral form is 14], 



*(ri, r 2 , ...r N ) = *(r 2 , r 2 , r 3 , r N )(l + r 12 /2) + r 12 ■ C(r 2 , r 3 , r N ), (5) 
where C(r 2 , r 3 , r^) is an unknown vector. The wave function functional also satisfies the 



electron-nucleus cusp condition which is 



V>(r, r 2 , ...rjv) = ^(0, r 2 , ...rjv).(l - Zr) + r ■ a(r 2 , ...r N ), (6) 
for a = 2. Here again a(r 2 , ...rjy) is also an unknown vector. 

In terms of the Hylleraas coordinates, the Hermitian single-particle operators noted 
above and the normalization operator may be expressed as W(s,t) where, respectively, 
W(s,t) = (s 2 +t 2 )/2, W(s,t) = j^, W(s,t) = + ^1], and W(s,t) = 1. In 

general, observables can be represented by single-particle operators expressed as W(s,t). 
The expectation of the operator W(s, t) which is 

fnr\x]w(3,t)*\x]dT 

{ 1 f**\x)*\x)dT ' [> 

can on substitution of the wave function functional ip[x] of Eq.(2) be written as 

J |$(«, (3; s, t)\ 2 [W(s, t)— <W>] [f(q; s, t, u) - 2f(q; s, t, u) + l]dr = 0. (8) 

Equivalently, Eq.(8) may be rewritten as 

POO 

/ e- 2as g(s)ds = 0, (9) 
Jo 

where 



g(s)= / duu dtcosh 2 (/3t)[W(s,t)- <W >](V - t 2 )[f 2 (q; s,t,u) - 2f(q; s,t,u) + 1}. 
Jo Jo 

(10) 

We now assume that the expectation (W) is known either through experiment or via some 
accurate calculation Q|. 

The next step is the constrained search over functions x(Q} s ) f° r which the expectation 
(W) of Eq.(7) is obtained. If the parameter a in Eq.(9) is fixed, then there exist many 



functions g(s) for which the expectation (W) can be obtained. This corresponds to a large 
subspace of wave function functionals (See Ref. 1). On the other hand, if the parameter a 
is variable, then the only way in which Eq.(9) can be satisfied is if 

g(s) = o, (li) 

This is equivalent to the constrained search of all wave function functionals over the 
subspace in which Eq.(9) is satisfied. 

Substitution of f(x',s,t,u) into Eq.(ll) leads to a quadratic equation for the function 

a(q, s) X (q; s) 2 + 2b(q, s) X (q; s) + c(q, s) = 0, (12) 

where 



ps pu 

a(q,s)= duu(l+u/2) 2 (l + qu) 2 e~ 2qu dtcosh 2 ([3t)(s 2 - t 2 ){W(s,t)- < W >], (13) 
Jo Jo 



b(q,s) = - du(l+u/2)(l+qu){e- 2qu (l+qu)-e- qu ] / dtcosh 2 ((3t)(s 2 -t 2 ) [W(s, t)- < W >], 
Jo Jo 

(14) 

c(q,s)= / du[e- 2qu (l + qu) 2 -2e- qu (l + qu) + l] dtcosh 2 (/3t)(s 2 -t 2 ){W(s,t)- < W >]. 
Jo Jo 

(15) 

Thus, in order to ensure that the wave function functional ip[x] leads to the exact expecta- 
tion value < W(s,t) >, one has to solve a quadratic equation for the determination of the 
functions x(<?; s). The subspace thus corresponds to two points. The two solutions Xi(<l', s ) 
and X2(<?; s) lead to two normalized wave functions ip[xi] and ip[x2\ each of which in turn 
give rise to the exact expectation < W(s,t) >. 

For the two normalized wave function functionals as determined above, the energy func- 
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tional in terms of Hylleraas coordinates which is 

ip*HiJ)dT (16) 



= Jo Jo Jo {U(S ^ + ^ +( ^ } ] 

-[AZsu- (s 2 -t 2 )]^ 2 }, (17) 
is then minimized with respect to the parameters a, /? and q. 



The above framework presented for the ground 1 X S state of the two electron system 
is general and also applicable to excited states. For example, if one were to consider the 
excited 2 3 S triplet state of the Helium atom, one could employ for the prefactor in Eq.(2) 
for the wave function functional ip[x) the expression $(a;s, £) = -y/|(~) e_Q ^- Note that 
in this simplest of choices used for explanatory purposes, screening effects are ignored. 
With such a choice, the procedure to determine the wave function functional ip[x] is the 
same as described above. In addition, this procedure could be employed in conjunction 
with the theorem of Theophilou a] according to which if ipi, if 2, , <Pmy> are orthonormal 
trial functions for the m lowest eigenstates of the Hamiltonian H, having exact eigenvalues 
Ei, E 2 , -Em,... , then Y1T= i(v 9 i|-^lv 9 «) — Y^Li-^i ■ m this way, a rigorous upper bound 
to the sum of the ground and excited states is achieved. With the ground state energy 
known, a rigorous upper bound to the excited state energy is then determined, while si- 
multaneously a physical constraint or sum rule is satisfied or an observable obtained exactly. 



The description of the constrained-search — variational method given in this section 
concerns the determination of wave function functionals that obtain the expectation value 
of arbitrary Hermitian single-particle operators exactly. The functions \ were assumed to 
depend only on the Hylleraas coordinate s, and as a consequence, a quadratic equation had 
to be solved for their determination. If the variational space is expanded, then one would 
have to solve an integral equation for the function x- 



The ideas of the constrained-search — variational method may also be applied to sum 
rules involving two-particle properties. For example, consider the pair-correlation density 
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g(rr') which is the conditional density at r' of all other electrons, given that one electron is 
at r, and which accounts for electron correlations due to the Pauli exclusion principle and 
Coulomb repulsion. The pair-correlation density for an N-electron system is defined as 

<7(rr') = <tf| - " r)|*>/p(r), (18) 

and satisfies the sum rule 

J g{vr')dr' = N - 1, (19) 

for each electron position r. However, in order to determine the wave function functional 
ip[x\ that satisfies this sum rule at each electron position, one must solve an integral 
equation for \. The details of the calculation of such a wave function functional are to be 
presented elsewhere jl5| . 



III. APPLICATION TO THE GROUND STATE OF THE HELIUM ATOM AND 
THE NEGATIVE ION OF ATOMIC HYDROGEN 

In this section we apply the constrained-search — variational method as described above 
to the ground state of the Helium atom and the negative ion of atomic Hydrogen. The 
constraint employed is that of normalization, and the prefactor is that of Eq. (4). We begin 
with a discussion of the wave function functional determined. 

Wave function junctionals 

The 3-parameter wave function functionals are determined by solution of the quadratic 
equation Eq.(12). This solution for the functions \ is analytical so that the wave function 
functionals ip[xi] and VDfca] too are analytical. We do not provide here the analytical 
expressions for xi(g,a,/5;s) and X2(g,a, P; s), but these functions are plotted in Fig.l. 
Observe that the two solutions for both He and H~ are distinctly different: one is positive 
and monotonically decreasing while the other is negative and monotonically increasing. 
Thus, although the two wave functions have the same structural form, and both satisfy 
the normalization constraint and the electron-electron cusp condition, they are very different. 
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FIG. 1: The functions Xl(?> a i A s) and X2(q,ct, (3; s) for H and He. 



The results as determined by these two wave functions for the ground state energy, 
and various single- and two-particle expectations are given in the subsections below. 
Comparisons are made with the results of the prefactor, Hartree-Fock (HF) theory, the 
3-parameter Caratzoulas-Knowles (CK), and 1078-parameter Pekeris wave functions. 
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Ground-state energy 

In Table I, we quote the values for the ground-state energy for H~ and He. The corre- 
sponding satisfaction of the virial theorem, and percent errors when compared to the values 
of Pekeris for He and those of the variational-perturbation results of Aashamar [16j for H 
are also given. Observe that the energies obtained by each wave function functional for H 
and He are an order of magnitude superior to that of the prefactor. For He, these results are 
on the average 0.06% from the Pekeris values. They are also an order of magnitude superior 
to both those of HF and CK. For H~ , both wave function functionals lead to results within 
0.1% of the Aashamar values, and to positive electron affinities as must be the case since 
the ion is stable. (In the HF approximation, one does not obtain the negative ion of atomic 
Hydrogen to be stable. The exact satisfaction of the virial theorem by HF theory, however, 
is a consequence of self-consistency.) The results clearly demonstrate that highly accurate 
ground state energies can be obtained by constructing few-parameter wave functions that 
are functionals. These energies are far superior to those determined by similar wave func- 
tions with the same number of parameters but ones that are not functionals. 

Single-particle expectations 

In this subsection we present the results of the expectations of the Hermitian single- 
particle operators W = J2i r i'y n = ~ 2,— 1,1,2, W = ^i^( r i)> an d W = X^<K r « — r )- 
We begin with the determination of the electron density p(r), which is the expectation of 
the operator W = ^ 5(r$ — r), and from which all the other single-particle expectations 
may be obtained. (Of course, these expectations may also be determined directly from the 
wave function functionals.) The density p(r) is also required for the determination of the 
nonlocal Coulomb hole charge distribution p c (rr') as explained in the following subsection. 

Now the wave function functionals are in terms of the Hylleraas coordinates (s, t, u) which 
involve the position of both the electrons or both their radial distances from the nucleus. 
The electron density p(r), on the other hand, depends only on the coordinates of one of the 
particles. Its determination from wave functions that are written in terms of the Hylleraas 
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TABLE I: Rigorous upper bounds to the ground state energy of H~~ and He in atomic units as 
obtained from the wave function functionals determined via the constrained-search — variational 
method, together with the values due to Hartree-Fock (HF) theory J], and the Caratzoulas- 
Knowles(CK) [lfl| . Pekeris f| and Aashamar wave functions. The satisfaction of the virial the- 
orem, and the percent errors compared to the values of Pekeris and Aashamar are also given. 



Ion or Atom 


Wave function 


Parameters 


Ground state energy 


0/ 

7o error 


—V/T 


ii- 




„, n fifil 9 R fl "}77Q7 


— U.OlOOU 


9 7"Z7A 


9 nnm 

Z.UUU1 




q = 0.274, a = 0.659, (3 = 0.308 


-0.52721 


0.1023 


1.9979 




q = 0.094, a = 0.659, (3 = 0.306 


-0.52712 


0.1194 


1.9982 


Aashamar 




-0.52775 




2.0000 


He 




a = 1.68585, (5 = 0.49732 


-2.87566 


0.9663 


2.0000 




q = 0.957, a = 1.662, (3 = 0.399 


-2.90158 


0.0736 


1.9975 




q = 0.242, a = 1.663, (3 = 0.399 


-2.90232 


0.0482 


1.9988 


HF 




-2.86168 


1.448 


2.0000 


CK 




-2.89007 


0.470 


1.9890 


Pekeris 




-2.90372 




2.0000 



coordinates is as follows. The electron density 

p(r) = J ij)*(^5(T i -T))il)dT = 2 J tp 2 (r r')dr', 

Using the symmetry of the two electronic system, we have 

r' 2 dr' / dcosO. 



With u = \/t 2 + r' 2 — 2rr'cos8, then, for fixed r and r we can rewrite Eq.(21) 



as 



co / r-r+r' 

dv' = 2n I —dr' / udu. 

r J\r-r'\ 



(20) 



(21) 



(22) 



On rewriting the wave function in terms of (r,r',u), and substituting Eq.(22) into Eq.(20) 
leads to 



f-oo / rr+r' 

p(r) = 2 I —dr' I uip 2 (r,r' ,u)du 

r J\r-r'\ 







Po(r) + Ap (r), 
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(23) 



where po( r ) is the density due to the prefactor (see the Appendix for the analytical expres- 
sion): 

p (r) = 2N 2 J e- 2as cosh 2 ((3t)dr', (24) 

and 

Apo(r) = 2N 2 J e~ 2as cosh 2 ((3t)(f 2 (x;s,t,u)-2f( X ;s,t,u))dr', 

(25) 

is the density due to the correlation term, which can be evaluated numerically. 
The electron density at the nucleus is 

p(0) = f pCEsinWdr 

J i 

= p (0) + Apo(O), (26) 
where po(0) is the prefactor contribution(see Appendix): 

p (r) = 2N 2 J e- 2ar cosh 2 ([3r)dr, (27) 
and the correlation contribution is 

Apo(r) = 2N 2 [ e- 2ar 'cosh 2 (Pr')(f(x; S ,t,u)-2f(x-,s,t,u))\ ri=r=0 , u= r 2 =r'dr'. (28) 



In Table II we quote the expectations of the operators W = 'Yj^^n = —2,-1,1,2, 
and W = ^j<5(i"i), for the ground state of the He atom as determined by the functionals 
ip[Xi] and ij)[x2\ together with those of Hartree-Fock theory, and the Caratzoulas-Knowles 
and Pekeris wave functions. The corresponding percent errors relative to the values of 
Pekeris are given in Table III. As expected (see Table III), the improvement over the 
prefactor values is significant. The results of the two wave function functionals and those 
of Hartree-Fock theory are essentially equivalent, indicating thereby that the corresponding 
densities are also essentially the same. The expectations of single-particle operators in 
Hartree-Fock theory are, of course, known to be correct to second order [gj]. Hence, both the 
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TABLE II: The expectation value of the operator W = Y2l=i r i"i n = — 2, —1,1,2 and 
W = Yli=i ^( r «) f° r the He atom employing the wave function functionals determined by 
the constraint-search — variational method, and by the Hartree-Fock theory(HF) 4], Caratzoulas- 
Knowles(CK)[lo|, and Pekeris |6j wave functions (WF). 



WF 


(*(n) + *(r 2 )) 


< (1/n + l/r 2 ) > 


< (1/rf + l/rl) > 


< ri + r 2 > 


< ( r i + r l) > 




3.5025 


3.3717 


11.930 


1.8758 


2.4757 


V'lxi] 


3.6295 


3.3750 


12.033 


1.8652 


2.4156 




3.6450 


3.3735 


12.048 


1.8639 


2.4112 


HF 


3.5964 


3.3746 


11.991 


1.8545 


2.3697 


CK 


3.3245 


3.3911 


11.714 


1.7848 


2.1292 


Pekeris 


3.62086 


3.3766 


12.035 


1.8589 


2.3870 



wave function functionals are accurate throughout space including the deep interior and 
far exterior of the atom. The comparison with the Caratzoulas-Knowles values (see Table 
III) is interesting for its implications. The wave function functional values are an order of 
magnitude superior. Of course, one does not expect the CK results to be accurate because 
these single-particle expectations are correct only to first order in the accuracy of the 
wave function. Thus, our results once again demonstrate, that wave function functionals 
determined by the constrained-search — variational method are superior to variationally 
determined wave functions that are not functionals. 



Structure of Coulomb holes 



We next consider the structure of the Coulomb hole charge distribution p c (rr ; ) as a 
function of the electron position r. The definition of this nonlocal or dynamic charge whose 
structure changes with electron position for nonuniform electron gas systems derives from 
that of the pair-correlation density g{vr') of Eq.(18) and from local effective potential energy 



theory 



The pair-density may be separated into its local and nonlocal components as 



g{rr') = p(r') + p xc (rr') 



(29) 



15 



TABLE III: The percentage errors of the results in Table II relative to the values of Pekeris. 



WF 


(<5(n) +5(r 2 )) 


< (1/n + l/ro) > 


< (1/r? + 1/ri?) > 


< 7"i + T2 > 


< (V? + Tn) > 




3.270 


0.145 


0.872 


0.909 


3.716 




0.237 


0.047 


0.017 


0.339 


1.198 




0.667 


0.092 


0.108 


0.269 


1.014 


HF 


0.676 


0.059 


0.366 


0.237 


0.725 


CK 


8.185 


0.429 


2.667 


3.986 


10.800 



where p xc (rr') is the Fermi-Coulomb hole charge. This dynamic charge distribution is the 
change in the pair density relative to the density that occurs as a consequence of the Pauli ex- 
clusion principle and Coulomb repulsion. It follows from Eq.(19) that its total charge is — 1. 
The definition of the Coulomb hole p xc (rr') derives in turn from that of the Fermi-Coulomb 
p xc (rr') and Fermi p x (rr') holes, the latter being defined through local effective potential en- 
ergy theory. In this theory, the interacting system as described by the Schrodinger equation 
is replaced by one of noninteracting Fermions with the same density. The corresponding 
wave function is a Slater determinant of single-particle spin orbitals, and one can then write 
down the resulting pair-correlation density p s (rr') of the model system as 

S s (rr')=p(r')+p,(rr'), (30) 

where p x (rr'), the Fermi hole, is the nonlocal component of this pair density, and is a 
consequence solely of the Pauli principle. The total charge of the Fermi hole is also —1. The 
Coulomb hole is then defined as the difference between the Fermi-Coulomb and Coulomb 
holes: 

Pc(rr') = p xc (rr') - p x (rr'), (31) 

and is thus representative solely of Coulomb correlations. The total charge of the Coulomb 
hole is 0. For two-electron systems in local effective potential theory 
then p x (rr') = — p(r')/2 independent of electron position r. 



the Fermi hole is 



In Figs. 2-4, we plot cross sections of the Coulomb hole p c (rr') for different electron 
positions r as obtained via the functional 1PIX2} together with the 'exact' Coulomb hole 
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FIG. 2: Cross-section through the Coulomb holes determined from the wave function functional 
V>[X2] for electron positions at (a)r = (a.u.), and (b) r = 0.566 (a.u.). The corresponding 'exact' 
Coulomb hole|3] cross sections are also plotted for comparison. The electron position is indicated 
by the arrow. 
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FIG. 3: The figure caption is the same as in Fig. 2 except that the cross sections plotted are for 
electron positions at (c)r = 0.8 (a.u.), and (d) r = 1.0 (a.u.). 
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FIG. 4: The figure caption is the same as in Fig. 2 except that the cross sections plotted are for 
electron positions at (c)r = 1.5 (a.u.), and (d) r = 5.0 (a.u.). 
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determined by Slamet and Sahni[7(. (The electron, indicated by the arrow, is on the z axis 
corresponding to 9 = 0°. The cross section through the Coulomb hole plotted corresponds to 
6' = 0° with respect to the electron-nucleus direction. The graph for r' < corresponds to 
the structure for & = ir and r' > 0.) The electron positions are at r = 0, 0.566, 0.8, 1.0, 1.5, 
and 5.0 (a.u.). It is evident from these figures that the Coulomb holes as determined 
from the functional ip[x2\ closely approximate the exact results for electron positions 
throughout space: in the interior, within the atom, near its surface and outside the 
atom, and in the far asymptotic region. Note the cusp representative of the electron- 
electron cusp condition at the electron position which is indicated by an arrow in the figures. 



Two-particle expectations 



As a consequence of the accuracy of the dynamic Coulomb holes obtained, we expect 
the results for the expectation of two-particle operators to also be accurate. In Table IV we 
quote the values for the expectations of the operators u 2 , u,l/u, l/u 2 , where u = |rj — tj\, 
together with the Hartree-Fock and Pekeris values. The corresponding percent errors 
compared to those of Pekeris are given in Table V. Once again, the results are an order 
of magnitude superior to those of the prefactor, and are accurate for both functionals, 
although those due to ijj\x?\ are consistently superior (see Table V). Of course, as expected, 
the Hartree-Fock theory results are not accurate. 

If one were able to write the expectation of arbitrary operators O as functionals of 



the density: (O) = (if)[p] 
Hohenberg-Kohn theorem 



O 



ip[p\) = 0[p], as is possible in principle according to the 
, then it is in the expectation of two-particle operators that 
the small differences between the Hartree-Fock theory density and those of the two wave 
function functionals would be exhibited. 



IV. CONCLUDING REMARKS 

The idea of expanding the space of variations in variational calculations by writing the 
wave function as a functional of functions is appealing not only because the functionals lead 
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TABLE IV: The expectation value of the operators u 2 , u, 1/u, 1/u 2 , where u = |rj — r,-| as obtained 
by the constrained-search — variational method together with those due to Hartree-Fock theory J] 
and the Pekeris [a] wave function. 



Wave function 


< u 2 > 


< u > 


< 1/u > 


< 1/u 2 > 


HF 


2.3694 


1.3621 


1.0258 


1.8421 




2.4757 


1.3957 


1.0000 


1.6998 




2.5325 


1.4271 


0.9388 


1.4300 




2.5236 


1.4241 


0.9434 


1.4532 


Pekeris 


2.5164 


1.4220 


0.9458 


1.4648 



TABLE V: The percentage errors of the results of Table IV relative to the values of Pekeris 



Wave function 


(u 2 ) 


(u) 


(1/u) 


(1/u 2 ) 


HF 


5.845 


4.217 


8.458 


25.764 


$ 


1.618 


1.853 


5.730 


16.045 




0.636 


0.355 


0.738 


2.373 


-to] 


0.286 


0.142 


0.254 


0.791 



to more accurate upper bounds for the energy with fewer parameters, but also because, as 
demonstrated in this work, they lead to wave functions that are accurate over all space. 
Thus, both single- and two-particle expectations are also determined accurately. Certainly, 
one could claim by comparison with the results of Hartree-Fock theory, but without rigorous 
proof, that single-particle expectations obtained thereby are correct to second order in 
the accuracy of the wave function. It is also evident that the accuracy of two-particle 
expectations lies somewhere between first and second order. In contrast, variationally 
determined wave functions that are not functionals are accurate only in those regions of 
space contributing to the energy. Thus, for such wave functions, it is the expectation value 
of only those single- and two-particle operators that appear in the Hamiltonian that are 
reasonably accurate. All other expectations are correct only to first order. 
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The results of the present work could be further improved as follows: by expanding the 
space of variations through the function x; by employing other more efficacious choices 
for the analytical form of the correlation factor and thus of the wave function functional; 
and by improving the prefactor. In our work so far, we have employed analytical forms 
for the prefactor. ( The results of our prefactor for the ground state energy of both 
H~ and He are superior to those of Hartree-Fock theory, see Table I.) Of course, one 
could employ the Hartree-Fock theory Slater determinant as the prefactor. Or one could 
employ a determinantal prefactor based on the orbitals generated within the local effective 
potential framework of Quantal density functional theory (Q-DFT). In principle, these 
orbitals generate the true electron density via a model system of noninteracting Fermions. 
The corresponding local potential within Q-DFT depends upon the wave functions of the 
interacting and noninteracting systems. Therefore, the corresponding orbitals generated 
are representative of electron correlations due to the Pauli exclusion principle, Coulomb 
repulsion, and the correlation contributions to the kinetic energy. 



Finally, we are presently investigating the use of wave function functionals in conjunction 
with Q-DFT for the many-electron case of N > 2. In these calculations, the antisymmetric 
determinantal correlated wave function functional employed is of the form 

V[x] = $ {0,,}n^(l - f(x; Ti, rj)). (32) 

Here is a Slater determinant that defines the state of the system and whose orbitals <pi 

are generated via the differential equation of Q-DFT, fix; r^, is a spinless correlation func- 
tional: f(x, r h rj ) = e-P 2r2 [l - X {R)(1 + r / 2 )]> where r = r 4 - r,, R = r 4 + Tj ,/3 = qp l/? \R), 
q is a variational parameter, and x{R) is determined by the constraint of the Coulomb 
hole sum rule for each electron position. This wave function functional satisfies the 
electron-electron cusp condition. In this instance an integral equation is solved to 
determine the function x(R)- Further, the products of the correlation functional are limited 
to lowest order since higher order products of these factors are less significant Q|. The 
highest occupied eigenvalue of Q-DFT differential equation corresponds in principle to 
the negative of the ionization potential^]. The region that contributes principally to 
this eigenvalue is the asymptotic classically forbidden region of the atom. In Q-DFT, 



22 



the asymptotic structure of the effective potential is due solely to Pauli correlations, and 
can be determined exactly. This is because the contributions to the potential due to 



Coulomb correlations and Correlation-Kinetic effects decay more rapidly than (— l/r)|llj. 
so that the potential in this region arises only from the Fermi hole charge which is defined 
through the Slater determinant of the orbitals. Thus, accurate ionization potentials cab be 
obtained via the use of correlated-determinantal wave function functionals in conjunction 
with Q-DFT. These are variational-self — consistent calculations that lead to upper bounds 
for the energy while simultaneously satisfying a nonlocal physical constraint. We are also 
currently investigating the construction of wave function functionals of the form employed 
in the present work, but with the satisfaction of constraints other than that of normalization. 



APPENDIX: 

We give the analytical expressions for the normalization constant, the energy, and various 
single- and two-particle expectation values as determined by the prefactor wave function 

$ = Ne- as cosh{pt). (A.l) 

Normalization 

/*00 [>S j*S 

dr$ 2 = 2n 2 N 2 / dse~ 2as / dtcosh 2 {f3t) I duu{s 2 - 1 2 ) 

2 2 -2a 6 + 3c* 4 /? 2 - 3a 2 /? 4 + /? 6 
= N ^ ( 2a6(/?-a)3(a + /?)3 } = (A ' 2) 

Ground-state energy 

E n = I <&*H§dT 



2 „ a(/? 2 -a 2 )(10a 4 -lla 2 /? 2 + 5/? 4 ) 
a — 2Za + 



8(-2a 6 + 3a 4 /? 2 - 3a 2 /? 4 + /? 6 
/? 4 (3a 4 -3a 2 /? 2 + /? 4 ) 



-2a 6 + 3a 4 /? 2 - 3a 2 /? 4 + /? 6 ) 



(A.3) 
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Expectation values 

Po(r ) = (f(ri _ r) + ,(„ _ r)) = Ae*[I + + (^)l. < A ^) 

*<0) = (iW + = JV 2 7r[Jj + + ^j)], (A.5) 



//»00 /*S /*S 

drs$ 2 = 2ir 2 N 2 / dsse~ 2as / dtcosh 2 {{3i) / duu(s 2 - t 2 ) 

Ar2 2 , 3(2a 8 - 4a 6 /? 2 + 6a 4 /? 4 - 4g 2 /? 6 + I 8 ) 
= " * ( 2a^-a) 4 (a + /3) 4 } ' ^ 



(— + —)= [ dr—^—® 2 = 2n 2 N 2 f dsAse~ 2as f dt cosh 2 ((3t) f duu 
n r 2 J s 2 - t 2 J Q J J t 

= 2a. (A.7) 



f dT 8< < S 2 = 2n 2 N 2 f°° dse' 2 ™ f ' dt cosh 2 ((3t)4(s 2 + t 2 ) 

J \ s ~ t ) J Jo 

- 4a V 2 + 6a 2 /? 4 - 

a 4 (/3-a) 3 (a + /?) 3 



2 2 (4a 6 - 4a 4 /3 2 + 6a 2 /? 4 - 2/3 6 ) 
^ ff 777^ 3771 , ma ■ t A -«) 



(lf + rf > = / <fr (^ > = J¥V fte- r a °°*w + *>»-«■)■ 

J 2 Jo Jo 2 

2 2 3(-2a 10 + 4a 8 /? 2 - 10a 6 /? 4 + 10a 4 /? 6 - 5a 2 /? 8 + /? 10 ) 
~ a 8 (/3-a) 5 (a + /5) 5 ' ( } 



rf™$ 2 = 2iV 2 7r 2 / dse~ 2as / cos/i 2 (/3t)(s 2 -t 2 ) / u 2 du 
Jo Jo Jt 

(70a 8 - 126a 6 /3 2 + 209a 4 /? 4 - 140a 2 /? 6 + 35/3 8 ) 



16a(/3 2 - a 2 ) (-2a 6 + 3a 4 /? 2 - 3a 2 /3 4 + /3 6 ) ' 



dm 2 ® 2 = 2iVV / dse- 2as / dt cosh 2 ((3t)(s 2 - 1 2 ) / M 3 d« 
Jo Jo Jt 

6(-2a 10 + 4a 8 /? 2 - 10a 6 /? 4 + 10a 4 /3 6 - 5a 2 /? 8 + /3 10 ) 



a 2 (/? 2 - a 2 ) 2 (-2a 6 + 3a 4 /? 2 - 3a 2 /? 4 + /? 6 ) 
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(A.10) 



(A.11) 



drU 2 
u 



poo ps ps 

2N 2 n 2 / dse- 2as / dtcosh 2 ((3t)(s 2 -t 2 ) / du 
Jo Jo Jt 

a{(3 2 - a 2 )(10a 4 - lla 2 (3 2 + 5(3 4 ) 



8(-2a 6 + 3a 4 /? 2 - 3a 2 (3 4 + (3 e ) 



(A.12) 



(i-) = [ dr^ 2 = 2N 2 ir 2 f°° dse~ 2as f dtcosh 2 ((3t)(s 2 - t 2 ) [" -du (A.13) 
n 2 J u Jo Jo Jt u 

Eq.(A.13) can be evaluated numerically. 
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